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Abstract. I review the basic ideas behind lattice QCD calculations that involve charm 
and bottom quarks. I report on the progress in getting the correct hyperfine splitting in 
charmonium from lattice QCD. Some of the basic technology behind numerical lattice 
QCD calculations is explained by studying some specific examples: computation of the 
charm quark mass, and the calculation of /s. 



1 Introduction 

The B factories at KEK and SLAC are producing a wealth of new data on the 
decays of the B meson. One of the main goals of the current heavy flavor program 
is to check the CKM matrix formalism by measuring the matrix elements with 
sufficient accuracy. To convert the experimental data into information about the 
quarks requires the accurate computation of hadronic matrix elements. The best 
(and some would say only) way of computing the required matrix elements is to 
use lattice QCD. 

In this paper I will explain the generic features of lattice QCD calculations 
that involve heavy quarks. Heavy quark lattice calculations share many common 
features to continuum calculations, such as matching to effective field theories. 
However, the more general formalism of lattice QCD allows a richer set of tools 
beyond just using perturbation theory. As most lattice QCD calculations share 
generic features, I will work through an example of computing the charm quark 
mass to show the important parts of the calculation. I will then discuss the 
calculation of the fg decay constant. I assume the reader is already familiar 
with the basic ideas of heavy quark effective field theory jl],|| . 

The latest results on heavy quark physics from the lattice are reported in the 
reviews at the annual lattice conference H|J|j| ■ The contents of the proceedings 
of the lattice conference have been put on hep-lat for the past couple of years P,[7 | 
The longer reviews by Kronfeld ||, Davies p|,|lO[ and Flynn and Sachrajda JlT | 
contain other perspectives on heavy quarks on the lattice. Gupta gives a general 
overview of lattice QCD [jL2| . There have been recent (political) developments 
to set up a working group on producing a "particle data table" for lattice QCD 
results U|. 
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2 A brief introduction to numerical lattice QCD 

Most lattice QCD calculations start from the calculation of the correlator Cij(t) 
defined in terms of the path integral as: 

aM = i / d u f d^o(t)iO(o)U- SF - Sa (i) 



where Sp is the action of the fermions and Sq is the action of pure gauge theory. 

The path integral is regulated by the introduction of a four dimensional 
space-time lattice. A typical lattice volume would be 24 3 48. The path integral 
is evaluated in Euclidean space for convergence. The fermion action is 

S F = i>Mtp (2) 

where M is called the fermion operator, a lattice approximation to the Dirac 
operator. The quadratic structure of the fermion action in equation || allows the 
integration over the fermion fields to be done explicitly. 

Cij{t) = - J dUO{t) t O{0))det(M)e~ Sa (3) 

The det(M) term controls the dynamics of the sea quarks. The 0(t) operator 
controls the valence content of the state. For example an operator {0(t)) for a 
B meson would be: 

0(t) = 5(t) 76 g(i) (4) 

where b and q are operators for the bottom and light quarks respectively. The 
operators in the path integral are Wick contracted to form a combination of 
quark propagators inside the path integral over the gauge fields. 

c(t) = ^J dUb(t) l5 q(t)q(0h 5 b(0)det(M)e- SG (5) 

The physical picture for the expression in |B| is a B meson created at time and 
propagating to time t where it is destroyed. Figure |l| shows the propagation of 
the light and heavy quark in the vacuum. 

The path integral expression for the correlator in equation [l] is calculated 
using Monte Carlo techniques on the computer. The ideas are sophisticated 
variants of the Monte Carlo method used to compute integrals. 

The algorithms, usually based on importance sampling, produce N samples 
of the gauge fields on the lattice. Each gauge field is a snapshot of the vacuum. 
The QCD vacuum is a complicated structure. There is a community of people 
who are trying to describe the QCD vacuum in terms of objects such as a liquid 
of instantons (for example |l4|j ). The correlator c(i) is a function of the bottom 
(M(U (i))^ 1 ) and light quark M(U(i))~ l propagators averaged over the samples 
of the gauge fields. 



1 N 

C W ~ N £ fmU(i))- b \M(U(i))- 1 ) (6) 
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The quark propagator is the inverse of the fermion operator. In perturbative 
calculations the quark propagator can be computed analytically from the fermion 
operator. In lattice QCD calculations the gauge fields have complicated space- 
time dependence so the quark propagator is computed numerically using variants 
of conjugate gradient algorithms. 

The physics is extracted from the correlators by fitting the correlator to a 
functional form such as [?]. 

c !J W^°e-"" i + c^ + ... (7) 

To visually judge the quality of the data, the correlators are often displayed as 
effective mass plots 

c(t) 

m*t f{t) = M^^y) ( 8 ) 

An example of an effective mass plot (using the data generated for |l5| ) is in 
figure |2[ 

The computationally expensive part of lattice QCD calculations is generating 
the samples of gauge fields. The most expensive part of a lattice calculation is 
incorporating the determinant in equation |[ The SESAM collaboration Jl6| 
estimated that the number of floating point operations (Nfi op ) needed for nf 
=2 full QCD calculations as: 

N flop *(-) 5 (—r s (9) 

a arripi 

A flop is a floating point operation such as a multiplication or addition. The 
value of Nfiop represents amount of calculation required on the computer and 
even more importantly the cost of the computer required. 

In some sense equation |9| (or some variant of it) is the most important equa- 
tion in numerical lattice QCD. To half the size of the pion mass used in the 
calculations requires essentially a computer that is seven times faster. Equa- 
tion ^| is not a hard physical limit. Improved algorithms or techniques may be 
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Table 1. Typical parameters in recent unquenched lattice QCD calculations. 



cheaper. In fact the "Asqtad" fermion action designed by the MILC collaboration 
is already computationally cheaper |l7j than the cost estimates in equation |^. 

The cost formula in equation ^ is for the generation of the gauge configura- 
tions. Once the gauge configurations have been generated, correlators for many 
different processes can be computed using some generalization of equation ^]. 
This class of calculation can be carried out on a farm of workstations. The lat- 
tice QCD community are starting to create publicly available source code ]l8[ 
and gauge configurations |Ig|pCj . 

Table [l] shows the parameters of some recent large scale unquenched calcu- 
lations. It is not considered necessary to do lattice calculations with physical 
light masses Mps/Mv = M^/Mp ~ 0.18. The aim is calculate with light enough 
quarks so that chiral perturbation theory can be used to make contact with the 
experiment. Sharpe pE| estimates that going as light as Mps/My ~ 0.3 may be 
necessary. 

The high computational cost of the fermion determinant led to development 
of quenched QCD, where the dynamics of the determinant is not included in 
equation ^, hence the dynamics of the sea quarks is omitted. Until recently the 
majority of lattice QCD calculations were done in quenched QCD. When the 
dynamics of the sea quarks are included I will call the calculation unquenched. 

Lattice QCD calculations produce results in units of the lattice spacing. One 
experimental number must be used to calculate the lattice spacing from: 

« = am fatt/ m expt (10) 

As the lattice spacing goes to zero any choice of m* xpt should produce the same 
lattice spacing - this is known scaling. Unfortunately, no calculations arc in this 
regime yet. The recent unquenched calculations by the MILC collaboration p6[ 
may be close. 

Popular choices to set the scale are the mass of the rho, mass splitting between 
the S and P wave mesons in charmonium, and a quantity defined from the lattice 
potential called tq. The quantity r$ is defined by tq p7f . 

rg^L=1.65 (11) 

Many potential |27j models predict r ~ 0.5 fm. There is no perfect way to 
compute the lattice spacing. Although it may seem a little strange to use r$ 
to calculate the lattice spacing, when it is not directly known from experiment, 
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there are problems with all methods to set the lattice spacing. For example, to 
set the scale from the mass of the rho meson requires a long extrapolation in 
light quark mass. Also it is not clear how to deal with the decay width of the 
rho in Euclidean space. 



2.1 Fermion actions for light quarks 

The lattice QCD formalism has the quark fields on the nodes of the lattice. 
The gauge fields are SU(3) matrices and lie between the nodes of the lattice. 
There are a variety of different ways of writing a lattice approximation to the 
Dirac operator on the lattice. Gupta jl2| reviews the problems and possibilities 
of fermion actions on the lattice. Discussions between lattice gauge theorists over 
which lattice fermion action is best must seem to outsiders to have the flavour 
of fanatical religious discussions, with the lattice "community" breaking into 
various sects, and accusations of "idolatry" being flung around. In the continuum 
limit all the fermion actions should produce the same results. This is clearly a 
good check on the results. 

As a starting point I will consider the Wilson fermion action. 

S T = E( k E^(^-1)^(^x+ m -? :c+ m(7m + 1)^(^x}+?xV'x) (12) 



The tree level relation between n and the quark mass m is n — 2 (4+ m ) • 

There is a lot of effort in the lattice gauge community on designing new 
fermion actions for light quarks (with masses lighter than the strange quark 
mass). There has been a long standing concern about fermion doubling on the 
lattice [[l2|. There are a number of pragmatic solutions to the doubling problem. 
For example the action in equation [l2] breaks chiral symmetry with an 0(a) 
term. Chiral symmetry will be restored as the continuum limit (a — > 0) is taken, 
but lack of chiral symmetry at finite lattice spacing causes problems, such as the 
difficulty of reaching light quark masses. 

Our understanding of chiral symmetry on the lattice has increased by the 



rediscovery of the Ginsparg- Wilson relation 28 



Mj 5 + 75 M = aMj 5 M (13) 

where M is the fermion operator in equation || at zero mass. 

Lattice fermion operators that obey the Ginsparg- Wilson relation (equa- 
tion 13) have a form of lattice chiral symmetry P9| . Explicit solutions, such 
as overlap-Dirac |5(i| or perfect actions |3l) , to equation |l3| are known. Actions 
that obey the Ginsparg- Wilson relation are increasingly being used for quenched 
QCD calculations 132]. This class of actions have not been used for heavy quark 
calculations (see ||33| for some speculations). Domain Wall actions, that can 
loosely be thought of as being approximate solutions to the Ginsparg- Wilson 
relation are being used in calculations [ p4p^ ] of the matrix elements for the 
e'/e. Unfortunately, solutions of the Ginsparg- Wilson relation are too expensive 
computationally to be used for unquenched calculations ]3q|. 
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Table 2. Mass of the strange quark from CP-PACS |22| in the MS scheme at a scale 
of 2 GeV. 



A more pragmatic development in the design of light fermion actions is the 
development of improved staggered fermion actions j37],|2l| . This class of action is 
being used for unquenched lattice QCD calculations with very light quarks (see 
table 0) by the MILC collaboration. The improved staggered quark formalism 
is quite ugly compared to actions that are solutions of the Ginsparg- Wilson 
relation. It is not understood why calculations using improved staggered quarks 
are much faster |l7[] than calculations using Wilson fermions [T^| . 

The largest systematic study of light hadron spectroscopy in quenched QCD 
has been carried out by the CP-PACS collaboration Q. CP-PACS controlled 
the systematic errors by using a 0.1 - 0.05 fm, m w /m„ ~ 0.75 - 0.4, and 
box sizes greater than 3 fm. A summary of CP-PACS's results is in figure || CP- 
PACS summarize their calculation of the light hadron spectrum in quenched 
QCD by the masses showed a deviation from experiment of less than 11%. 

Although an agreement between experiment and the results of quenched QCD 
at the 11% level might seem impressive, many heavy-light matrix elements need 
to computed to an accuracy of under 5% to have an impact on tests of the CKM 
matrix. 

There are now indications of the effects of the sea quarks. For example, 
the CP-PACS collaboration j22) have used nj — 2 lattice QCD to calculate the 
strange quark mass. CP-PACS's results are in table The results show a sizable 
reduction in the mass of the strange quark between quenched and two flavour 
QCD. See the review by Lubicz |39| for a discussion of the results of CP-PACS 
in comparison to those from other groups. The results from CP-PACS for the 
strange quark mass need to be checked by other calculations with lighter quarks 
and finer lattice spacings, 



3 The different ways of treating heavy quarks on the 
lattice 

In principle all the above formalism can be used to do calculations that include 
charm and bottom quarks. Unfortunately, in practice there is a restriction that 
the quark mass should be much less than the lattice spacing. 

aM Q « 1 (14) 

As the lattice spacings accessible to current calculations are - ~ 2 Gev (see 
table [l]) and 4 GeV Q] for unquenched and quenched QCD respectively. Hence 



Focusing of a Parallel Beam 7 



calculations using traditional techniques will just about work for the charm mass 
(m c ~1.3 GeV), but will not work for the bottom quark (m& ~5 GeV). It is 
computationally expensive (see equation ||) to reduce the lattice spacing, so that 
a b quark will be resolved by the lattice. There are a variety of special techniques 
for including the bottom quark in lattice QCD calculations, all of them arc based 
on heavy quark effective field theory. 

I do not discuss the method of using the potentials measured in lattice 
QCD with Schrodinger's equation to compute the mass spectrum of heavy-heavy 
mesons. This subject is reviewed by Bali EJ. The potential based approach is 
not applicable to computing matrix elements, so is not useful for checks of the 
CKM matrix. 

Various subsets of the lattice QCD community have strong opinions on the 
right (and wrong) approach to including heavy quarks in lattice calculations. 
Obtaining consistent results from calculations that use different heavy quark 
actions is a good check on the systematic errors. In the next sections, I describe 
some of the more popular techniques used for heavy quarks on the lattice. I 
report on the results from this class of methods in section 3.5. 



3.1 The improvement view 

There are concerns that the results from lattice calculations have large lattice 
spacing errors because ciMq ~ 1 even for charm quarks. The Wilson action in 
equation [l^ has 0(a) lattice spacing errors. If the 0(a) term is removed then 
perhaps larger quark masses could be used in the lattice QCD calculations. 

A standard technique from numerical analysis is to use derivatives that are 
closer approximations to the continuum derivatives. For example the lattice 
derivative in equation [l6| should be more accurate with a larger lattice spac- 
ing than derivative in equation [Po] , 

f(x + a) - f(x - a) df 2 

2a = T X + ° {a) (15) 

4 f(x + a)-f(x -a) _ f(x + 2a)-f(x-2a) = ±, q( as (w) 
3 1 2a 16a ' dx {G ' { 1 

However in a quantum field theory there are additional complications, such as 

operators mixing under renormalization. 



There is a formalism due to Symanzik 1 42 43fl called improvement where new 
terms are added to the lattice action that cancel 0{a) terms in a way that is 
consistent with quantum field theory. The "simplest" improvement []44J1 to the 



Wilson action is to add the clover term [Tt] to remove tree level lattice spacing 
errors: 

gclover = S W + ^(^<7^i^*) (17) 

X 

where is the lattice field strength tensor. 
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If the csw coefficient is computed in perturbation theory is used then the 
errors are 0(a<? 4 ). The ALPHA collaboration [Q have computed csw to all 
orders in g 2 using a numerical technique. The result for csw from ALPHA is: 

1 - 0.656 5 2 - 0.152 5 4 - 0.054 5 6 
CSW = 1 - 0.922^ ^ 

for < g < 1, where g is the coupling. The estimate of csw, by ALPHA 
collaboration, agrees with the one loop perturbation theory for g < 1/2. 

Some groups tried to use the results from lattice QCD calculations with 
quark masses around charm with the scaling laws from heavy quark effective 
field theory to compute matrix elements for the b quark. An example tried by 
the UKQCD collaboration was to extrapolate the Jb decay constant from 
masses around charm, where the clover action can be legitimately used, to the 
bottom mass, using a functional form based on HQET [El. 



The review by Bernard j4j describes the potential problems with this approach. 
An extrapolation in mass from 2 GeV to around the bottom quark mass at 
5 GeV is problematic. Note that UKQCD |Q did address some of Bernard's 
criticism || 

Computers are fast enough to directly include quark masses close to the 
bottom mass in quenched calculations. However, this approach will not work for 
unqucnchcd calculations for some time. 



3.2 The static limit of QCD 

It would clearly be better to interpolate in the heavy quark mass rather than use 
extrapolations. The static theory of QCD can be used to compute the properties 
mesons with light ante-quark and static (infinitely heavy) quarks. A combined 
analysis of data from static-light and heavy-light calculations can be used to 
interpolate to the b quark mass. This was the approach taken by the MILC 
collaboration |l7| in one of the largest calculations of the fs decay constant. 
The lattice static theory of Eichten and Hill [ [48[ 

S s tauc = ia 3 b\x)(b(x) - U a (x ~ t)b{x - t)) (20) 

X 

has been used for B meson physics. 

One of the reasons that static quarks have not been included in many calcu- 
lations is that it can be difficult to extract masses and amplitudes using equa- 
tion |?j because the signal to noise ratio is poor. However there are numerical 
techniques |49|| that are better, but not in wide spread use in matrix element 
determinations . 

To extract matrix elements from static-light calculations requires the static- 
light operators to be matched to QCD. The ALPHA collaboration have started |5(| 
a program to compute the matching and renormalization factors numerically. 
The one loop matching factors are available. 
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3.3 Nonrelativistic QCD 

It would be better to actually do lattice calculations at the physical bottom 
or charm quark masses, rather than extrapolate or interpolate to the physical 
points. The formalism called nonrelativistic QCD(NRQCD) allows this [5l|j5^| . 
NRQCD is an effective field approximation to QCD for heavy quarks. The op- 
erators in the NRQCD Lagrangian are ordered by the velocity v. 
A low order Lagrangian for NRQCD is 

3 

C = Q(At ~ E " Cnr TaT~ + ( 21 ) 

i r^ J ZMqcl ZMqci 

where are the covariant derivatives on the lattice and B is the chrmomagnetic 
field. 

The NRQCD formalism works both for both heavy-light (B) and heavy- 
heavy systems (T). Estimates from potential models |53| suggest that the v 2 ~ 



0.3 in charmonium. In section 3J5, I review the evidence that 
shows that the NRQCD is not as convergent in charmonium as the naive power 
counting arguments suggest. 

The main theoretical disadvantage of NRQCD is that the continuum limit 
can not be taken because of the ^ terms in the Lagrangian. In practice im- 
provement techniques can be used. 

The coefficients, such as c^r in equation |2l|, in the Lagrangian are fixed by 
matching to QCD. The matching calculations involve lattice perturbation theory 
that is harder than continuum perturbation theory because the Feynman rules 
are more complicated |54| . 

A physically motivated (but not rigorous) way of improving the convergence 
of lattice perturbation theory is to use tadpole improvement [p3| . Tadpole per- 
turbation theory can be used to produce "reasonable" tree level estimates for 
coefficients in the Lagrangian. 

To do perturbation theory the gauge links are expanded in terms of the gauge 
potential: 

U^x) = e'^'W -> 1 + iagA^x) + ... (22) 

Equation p2| suggests that the (U) ~ 1, however this is not seen in numerical 
lattice calculations. Also the complicated vacuum structure of QCD would make 
(U) ~ 1 unlikely. Lepage and Mackenzie j53) suggest: 

U^x) = e'^'W - u (l + iagA^x)) (23) 

where u is the "mean gauge link" . Unfortunately there is no unique way of 
defining the mean gauge link. The expectation value of the mean link is zero 
because it is not gauge invariant. Some estimates are based on taking the quartic 
root of the plaquette 

u .p = (^ReTriUpiag)) 1 ^ (24) 
or computing the mean link in Landau gauge. 
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There are projects p5j under way that attempt to estimate the coefficients 
such as cat 7? to order a 2 . The basic idea [|56| is to try to use weak coupling numer- 
ical lattice QCD calculations to obtain information on perturbative quantities. 

A similar approach to NRQCD is taken by the Fermilab group |37],[38| , except 
that they match to a relativistic fermion action, essentially the clover action with 
mass dependent coefficients. 



3.4 Anisotropic lattices. 

A technique that has been used for heavy quarks p9| , p0pl| | is to use a lattice 
spacing that is smaller in the time direction than in the space direction to cir- 
cumvent the restriction arriQ << 1. A finer lattice spacing is used in the time 
direction such that atm q < 1 but a larger lattice spacing a s is used the spatial 
direction to keep the cost down and stop any problems with finite size effects. 
This approach assumes that the discretization error is only weakly dependent 
on a s m q . 

The anisotropic clover operator |5l] is 

At = mo + voWojo + ^^Wai + l(w ^aoiFbi(x) + ^ V ^F io {x)) (25) 
to & — to . 



The clover action (equation [12] plus 17) is reproduced by the conditions: £o = 1, 
v = Vq, and w = Wq. 

The parameters: wq, u>i, vq, and v need to be correctly chosen. For example 
Klassen |5i| proposed to tune v, by computing the pseudoscalar meson mass at 
nonzero momentum, and to choose the value of v that gave c(p) = 1. 

E(p) 2 = E(Q) 2 + c(p) 2 p 2 (26) 

Klassen's J59| original motivation for using this class of action was that is was po- 
tentially easier to tune the unknown parameters using the techniques developed 
by the ALPHA collaboration |l5j than for the Fermilab heavy quark action |5TJ . 
The pure gauge action is also modified j6lj 

s 9 = P(r E [i-PsAx)] + ^YM- p ss{^)]) (27) 

to — , 

x,s>s' x,s 

where P ss are purely spatial plaquettes and P ss i are plaquettes in space and 
time. The renormalized anisotropy £o = a s /a t ( ratio of lattice spacings in time 
and space) can be measured by comparing the lattice potential in space and 
time ||. 

A practical problem in lattice calculations is that the signal in equation ^ is 
lost in the noise for large times. The smaller lattice spacing in the time direction 
from anisotropic lattices means that the region in lattice units, where there is a 
signal, is longer, thus it is easier to fit equation |^ to the data. Collins at al. Q 
used this feature of anisotropic lattices to get improved signals for form factors. 
Although the fit region in lattice units is longer, the actual fit region in physical 
units may be smaller, this may cause problems. 
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3.5 The hyperfine splitting in charmonium 

It is obviously important to test the methods used to solve lattice QCD by 
comparing the results against experiment. This validation procedure ensures 
that the various errors in the calculations are under control. Figure || shows the 
charmonium spectrum from lattice QCD calculations by the CP-PACS collabo- 
ration pH . The overall agreement with experiment is quite good. 

A particularly good test of lattice QCD techniques is to compute the mass 
splitting between the J/i/j and rj c . In the section I will use Amu to denote the 
mass difference between J/ip and r] c . The masses of these two meson can usually 
be computed with the smallest statistical error bars. Also, as these masses are 
independent of light valence quarks, this splitting does not depend on a chiral 
extrapolation of the valence quarks. Hein et al. |64| discuss the various systematic 
errors in the mass splittings between other heavy hadrons. 

I will start by discussing the results from quenched QCD. I will use Aran 
as the mass splitting between the J/ip and r\ c . The experimental value for the 
mass splitting between the J/tp and r\ c is 116 MeV. It was therefore disappoint- 
ing that some of the lattice QCD calculations in the early 90's were: Atuh = 
28(2) MeV and 52(4) MeV from the Wilson action and tree level clover action 
respectively |35| (a -1 = 2.73 GeV from the string tension). 

Using the clover action with the a value of the clover coefficient csw (see 
equation |l7| ) motivated by tadpole perturbation theory, the Fermilab group pc[ 
obtained Am H = 93 ± 10 MeV. 

The hyperfine splitting is sensitive csw &t nonzero lattice spacing, but the 
hyperfine splitting should be independent of the csw as the continuum limit is 
taken, because the clover term is an irrelevant operator. 

Recently the QCD-TARO collaboration |67j have studied the charmonium 
spectrum using the clover action at a smaller lattice spacing (a -1 ~ 5 GeV) 
than previously used. Using the clover action QCD-TARO collaboration p7| 
obtained Am H = 99 ± 7 MeV, AM = 87 ± 2 MeV, for quadratic and linear 
extrapolations in the lattice spacing to the continuum. 

The NRQCD collaboration calculated Am H to be 96(2) MeV fjf] using an 
NRQCD action that included terms of 0(m c v 4 ). The lattice spacing was a -1 
= 1.23(4) GeV. Using power accounting arguments the size of the next order 
in the NRQCD expansion was estimated to be 30 to 40 MeV. Unfortunately, 
when Trottier |i9) included the 0(m c v 6 ) relativistic corrections, he obtained 
AmH = 55(5) MeV. The hyperfine splitting is sensitive to the cnr coefficient 
(see equation ^j]). Further work by Shakespeare and Trottier j70| showed that 
the hyperfine splitting was sensitive to the tadpole prescription used to estimate 
cnr, so the final word on the utility of NRQCD for charm quarks may have to 
wait for cnr to be computed beyond one loop. The caveat is that the current 
estimates for cnr seem to produce good agreement with the hyperfine splittings 
in the baryon sector fn§ . 

Although NRQCD is clearly not the technique to use to compute the hy- 
perfine splitting in charmonium, NRQCD may be valid for hadrons with charm 
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quarks, such as the mass splittings between the S-wave states and the speculated 
1"+ state Q and D mesons @. 

There was a preliminary attempt ]73| to use an action motivated by renor- 
malisation group arguments (perfect action Q) to compute the charmonium 
spectrum. This class of action should produce results with reduced lattice spac- 
ing dependence. Unfortunately the action used in [j73| did not produce a result 
for the hyperfme splitting closer to experiment than any other approach. 

Klassen ]5j| first used anisotropic lattices to study the charmonium system. 
Using spatial lattice spacings in the range 0.17 to 0.3 fm (the scale was set 
using tq = 0.5/m) and two anisotropies 2 and 3, Klassen obtained a hyperfme 
splitting of just over 90 MeV in the charmonium system. Chen |}0| obtained 
Am H = 71.8(2.0) MeV with anisotropy £ = a s /a t =2. 

The definitive study of the anisotropic lattice technique for charmonium spec- 
troscopy was carried out by CP-PACS |3l) . They fixed the anisotropy at 3 and 
used spatial lattice spacings between 0.07 and 2 fm (finer then both Chen |Q 
and Klassen The results from CP-PACS |l) were: Am H = = 73(4) MeV 

using ro to set the scale and 85(8) MeV using the P-S splitting for the lat- 
tice spacing. I have combined the different errors using quadrature. CP-PACS 
concluded that a s m q < 1 is still required for a reliable continuum extrapolation. 

A qualitative explanation for the low value of the hyperfme splitting in char- 
monium from quenched QCD was given by El-khadra |6q| using potential model 
ideas. In El-khadra's model the Richardson potential p5[ is used 

V(q 2 ) = G F ^j 1 (28) 

/V <r 1 °g( 1 + 1 A ) 

with 

ffi = 11 - 2n//3 (29) 

to solve for the wave function of the charm quark. The wave function depends 
on the number of flavours n / . El-khadra obtained the result 

¥W) = °' 86 (30) 

In this model the hyperfme splitting is related to the wave function and coupling 

(a s ) as 

Am H ~ | W(0) | 2 (31) 

mi 

Including the suppression of the coupling in the quenched theory, El-khadra 
estimated 

^quenched _ 7QMe y (33) 

There have been a number of unquenched lattice QCD calculations |24 76[23| , |77f 
that have seen evidence for the n f dependence of the heavy quark potential at 
small distances. The MILC collaboration |77j have systematically studied the 
wave functions from the measured heavy quark potential from quenched and 
unquenched calculations. 
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There has not been much work on the charmonium spectrum from un- 
quenched lattice QCD calculations. El-Khadra et al. (78) did look at the charmo- 
nium spectrum on (unimproved) staggered gauge configurations from the MILC 
collaboration. No significant increase in the hyperfine splitting was reported. 
The m^/nfip was 0.6 and the lattice spacing was a" 1 ~ 0.99(4) GeV. 

Stewart and Koniuk [7S[ | studied the charmonium spectrum using NRQCD 
on unquenched (unimproved ) staggered gauge configurations {m n /m p ~ 0.45 
and a ~ 0.16 fm). Any signal for the effect of unquenching was hidden beneath 
the other systematic uncertainties in using NRQCD for charmonium. 

Although the potential model argument of El-Khadra for the effect of quench- 
ing on the hyperfine splitting gives some insight, it does not explain the sea quark 
mass dependence of the splitting. Grinstein and Rothstein |$0| have developed a 
formalism based on Chiral Lagrangian for the dependence of quarkonium mass 
splittings between IP-IS and 2S-1S on the sea quark mass. Up to chiral logs 
they predict for the splitting 5m 

5m ~ A + Bml (33) 

where A and B are unknown parameters and m n is the mass of the pion made 
out of light sea quarks. 

In my opinion a decade's worth of lattice QCD calculations of the J/ip -rj c 
mass splitting can be summarized as waiting for unquenched gauge configu- 
rations with light sea quark masses. Recently, progress has been made in the 
Upsilon system using the gauge configurations generated by the MILC collabo- 
ration. The preliminary work by Gray et al [jSl| found that the correct ratio was 
produced for the (P-S)/(2S-1S) mass splittings in Upsilon. It will be interesting 
to see the charmonium spectrum on these lattices, particularly if relations such 
as [53| can be tested and used. 

There is another possible reason that the hyperfine mass splitting between 
the J/ip and ry c is smaller than experiment in current simulations. All lattice 
calculations have computed the non-singlet correlator (see figure ^). However, 
charmonium interpolating operators are actually singlet, so the Wick contrac- 
tions contain bubble diagrams (see figure ||). The bubble diagrams are OZI sup- 
pressed so should be small. However, this argument will fail if there is additional 
nonperturbative physics. For light mesons p2), it has been found that the ef- 
fects of the bubbles can be large for the pseudoscalar and scalar mesons where 
the additional physics is the anomaly and the ++ glueball, but not for other 
channels. 

Morningstar and Peardon [[S3| have computed the glueball spectrum in quenched 
QCD. They obtained masses of 2590(40)(130) MeV and 3640(60)(180) MeV for 
the ground and first excited states of the h glueball respectively. Morningstar 
and Peardon computed the mass of the 1~~ glueball to be 3850(50)(190) MeV. 
So it is not inconceivable that the rj c mass (2980 MeV) is effected more by glue- 
ball states than the J/tp state. The above comments are speculations and can be 
checked by explicit lattice calculations. As the effect of the bubble diagrams is 
almost certainly less than 50 MeV, hence this will be a very hard mass splitting 
to estimate. 
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Grou 


E_ 


Method 


Afo* - Mc MeV 


Boyle | 
Boyle [ 
Hein et a 
UKQCD 


M 

~_ 

~m 

m 


clover 
/3=6.0 tadpole clover 
] NRQCD (3 = 5.7 
NP clover j3 = 6.2 


124(5)(15) 
106(8) 
103™ (3) (6) (5) 
130±«t» 


PDG [S6J 


Experiment 


140.64(10) 



Table 3. Collection of hyperfine splittings between the D and D* . 



The mass spectrum of heavy-light mesons introduces the additional compli- 
cation of the light valence quark. Lattice QCD calculations can be done with 
quark masses around charm, but for computational reasons the light quarks have 
masses that are typically greater than half the strange quark mass. The lattice 
data is extrapolated in the light quark mass to the physical points using a fit 
model based on chiral symmetry. 

In table || I have collected some results for the D* — D mass splitting. 
Currently there is a lot of effort in the lattice gauge theory community to study 
the chiral extrapolations of quantities with the light quark mass. The Adelaide 
group |87],^8| have developed various phenomenological forms for the light quark 
mass dependence of hadron masses loosely motivated by effective field theories. 
The fit models have had some empirical success with extrapolating the masses of 
the rho and nucleon |8?],^8| , with the caveat that lattice spacing errors were not 
taken into account. Similar techniques were applied to the B* — B and D* — D 
mass splittings by Guo and Thomas J8{|. No improvement with the agreement 
between experiment and the lattice data was seen. 



4 Case study: calculating the charm mass from lattice 
QCD 

The general steps involved in many lattice QCD calculations are fairly similar. To 
explain the component parts of a lattice QCD calculation, I will explain the use 
of lattice QCD data to extract the charm quark mass from experimental data. 
Full details of the calculation of the charm quark mass can be found in . I 

will not discuss the approach |j92| , j93| to computing the charm quark mass based 
on Fermilab's heavy quark action J57]]) and the pole mass. There is a useful 
review by El-Khadra and Luke |94[] on computing the mass of bottom quark. 

The error for the mass of the charm quark quoted in the particle data table is 
8% - an unbelievably large error for a basic parameter of the standard model that 
was discovered in 1974. The experimental mass of the D meson is 1869.3 ± 0.5, 
hence the error on the charm quark mass is predominantly from theory. There are 
many places in particle physics where a more accurate value of the charm mass 
would be useful. Some models of quark matrices predict relationships between 
quark masses and CKM matrix elements. For example Fritzsch and Xing p5[ 
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Kh 




aM P 


aM v 


0.1200 


0.1346 


0.841+} 


0.871+2 


0.1200 


0.1351 


0.823+1 


0.856+2 


0.1233 


0.1346 


0.739+} 


0.775±ij 



Table 4. Heavy-light meson mass from UKQCD J46J as a function of k value. Mps 
and My are the pseudoscalar and vector meson masses 



predict 

I Kb | = fm^ .v 
\Vd>\ V m c [ ' 

To test such relations we need to accurately determine all the component pa- 
rameters of the standard model. 

The starting point, I shall take is the masses of the heavy-light mesons as 
a function of lattice parameters. The hadron masses come from a fit of the 
correlator in equation For example, table ^ contains the masses of a heavy- 
light meson in lattice units j3 = 6.2 from UKQCD [H. The k value is defined in 
the action in equation [l2| After this point no supercomputers are required, just 
a nonlinear x 2 fitting program, physical insight and theoretical physics. 

To start the journey from the lattice to the real world, the lattice parameters 
need to be converted to more physical parameters. The first job is to convert from 
the kappa value into the quark mass. There are a number of different expressions 
for the quark mass in terms of the lattice parameters. 

One definition of the quark mass is based on the vector ward identity. 

1,1 1 N 

mV= 2 ( K-—? (35) 

fx nj crit 

The my quark mass suffers from an additive renormalisation for Wilson like 
fermions. The value of — - — is obtained by the value of k that gives zero pion 
mass. To remove 0(a) corrections the improvement formalism requires that 

m = m(l + b rn ma) (36) 

The perturbative value of b m is — ~ — 0.096g 2 |96[ . 

The quark mass can also be defined in terms of the PCAC relation. 

(a 4 At(t)P(o)) 

= wwm ( } 

There are also O(a) corrections to equation |37], see |91| for details. 

In principle the masses my and vhaw should agree, however at finite lattice 
spacing, where the calculations are actually done they disagree. For quark masses 
below strange, it has been shown that the two definitions agree as the lattice 
spacing is taken to zero (see |9j for a review). 

The masses must be converted from lattice units into MeV. As explained in 
section ^, one quantity must be sacrificed to find the lattice spacing. For example 
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at /3 = 6.2 in quenched QCD, UKQCD find a" 1 = 2.66±£, 2.91;}, 2.54+f, GeV 
from /„■, ro, and m p respectively. The spread in different choices should reduce as 
the continuum limit is taken in an unqucnchcd lattice QCD calculation. If there 
arc different choices of lattice spacing, this is usually included in the systematic 
error. 

Now we have a table of data heavy-light meson masses in GeV versus the 
lattice quark masses in GeV. The meson masses must be interpolated and ex- 
trapolated to the physical meson masses. The theory behind the extrapolations 
is an effective Lagrangian for mesons with heavy quark and chiral symmetry Q . 

The value of k corresponding to the strange quark mass is usually determined 
from light quark spectroscopy by interpolating to the mass of kaon or phi meson. 
Becirevic et al. investigated using three different fit models to extrapolate 
the meson masses in the heavy quark mass. 

M 1 HL (m Q )=A 1 +B 1 (m Q ) + C 1 m 2 Q (38) 

M 2 HL (m Q )=A 2 + B 2 (-^) + C 2 (-± T ) (39) 
iriQ nig 

M 3 HL (m Q ) = A 3 + B 3 (—) + C 3 m Q (40) 

The 1/Mq terms are motivated by heavy quark symmetry [^| . Kronfeld and 
Simone have used the fit model in equation ^ to estimate Ai and \ 2 pa- 
rameters of HQET. 



4.1 Quark mass renormalization factors 

In the last section I showed how to find the charm quark mass in the lattice 
scheme. However, normally quark masses are used in application in the M S 
scheme. Also, a consistent scheme is also required to compare the results from 
different calculations. 

The quark mass in the lattice scheme (mz,(a)) is matched to the MS scheme 

m Ms(^) = Z m (a^)m L (a) (41) 

where Z m [a[i) is the matching factor. 

The matching factor has been computed in perturbation theory to one loop 
order. 

Z m (a,j) = 1-^(8 ha(H - Cm) (42) 

The value of Cm depends on the fermion action. For the clover action, Cm — 
25.8 |§. 

As usual with one loop calculations, there is an ambiguity as to what scale 
(fi) to evaluate the matching at The "best guess scale" for the /i (called q-k) 
can in principle be computed using the formalism described by Lepage and 
Mackenzie |33| . Most people include the effect of varying fi in some range from 
1/a to 7T I a in the systematic errors. 
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The accuracy of the quark mass determination would improve if the matching 
could be done to higher order than one loop. The Feynman rules on the lattice 
are more complicated than in the continuum, hence calculations beyond one loop 
are very hard. Some groups are starting to try to automate lattice perturbation 



theory |100| 



The general framework of lattice QCD allows other approaches to computing 
matching factors without using Feynman diagrams on the lattice. Sint ]101[ and 
Sommer review some of the ways that matching factors are computed on 
the lattice. For example the a 3 term of the residual mass (important for the 
extraction of the bottom quark mass) of s tatic theory in quenched QCD was 



computed using a numerical technique 1 102 



A general technique [ 103 1 for matching between the lattice and MS schemes 
has been developed by the Rome and Southampton groups. The basic idea is to 
use the quark propagator calculated in lattice QCD to do the lattice part of the 
matching. The gauge has to be fixed in this approach. Usually Landau gauge is 
chosen. 

4.2 Evolving the quark mass to a reference scale 

The matching procedure produces the charm mass at the matching scale. To 
compare different mass determinations, the quark mass has to be evolved to a 
standard reference scale, essentially the same as that used by the particle data 
table. The reference scale chosen for the charm quark is the charm quark mass 
itself. 

The running quark mass equation is used to evolve the quark mass to the 
standard reference scale of the charm quark mass. 

M 2 i^(M) = m»'(Mb£'(a?') = -£ 7 ^(^)^ (43) 
The coefficients are known to four loop order. The required equations are 



conveniently packaged in the RunDcc Mathematica package [104]. The ALPHA 



collaboration have a method to do the evolution of the quark mass numeri- 



cally [105,59. The method is starting to be used for unquenched QCD [104 



4.3 Comparison of the results 

I have outlined the basic ideas behind a lattice QCD calculation of the charm 
quark mass. In table || I collect the state-of-the art results for the charm quark 
mass from (quenched) lattice QCD. These should be compared with the result 
quoted in the particle data table of m£* s (m c ) between 1.0 to 1.4 GeV. 

5 The f B decay constant 



A crucial quantity for tests of the CKM matrix formalism is equation 44 



Am s V ts 1 2 m B s c2 



Am d V t d m Bd 



e (44) 



18 Craig McNeile 



Group 


mf s {m c )GeV 


Becirevic et al. [M 
Rolf and Sint |lf 
Juge H 
Kronfeldgl 


1.26(4)(12) 
1.301(34) 
1.27(5) 
1.33(8) 



Table 5. The charm quark mass from quenched lattice QCD. 



N f 


Decay constant 





f B 


= 173(23) MeV 


2 


f B 


= 198(30) MeV 





fr>s 


= 230(14) MeV 


2 


fDs 


= 250(30) MeV 



Table 6. Summary || of the results for heavy-light decay constants from quenched 
Nf =0 and (Nf = 2) unquenched lattice QCD. 



The value of Am^ has been measured experimentally, while Am s is expected 



to be measured at run II of the Tevatron [106]. There are already useful exper- 
imental limits on Am s . The hard part is extracting the ratios of QCD matrix 
elements in 



1 4.",-. 



The quantity £ can not be extracted from experiment and is non-perturbative. 

The Jb decay is the QCD matrix element for the semi-leptonic decay of the 
B meson. It is analogous to the pion decay constant. It is claimed that Jb will 
never be measured experimentally, hence it must be computed from QCD. The 
B (bag) factors are also QCD matrix elements that have been computed from 
lattice QCD. 

The computation of Jb shares many features to the calculation of the charm 
quark mass. The same data from supercomputers could be used for both the fs 
and charm quark mass calculation, fs is extracted from the matrix element 

(0 | 4) | Qq,P = 0) = -if Qg M Qq (46) 

This matrix element is simply related to the amplitudes (cy) in equation 0. The 
main additional complication over the charm mass calculation is the extrapola- 
tion of the decay constant to the bottom mass. 

Sinead Ryan reviewed the latest results for the fs decay constant at the 
lattice 2001 conference ||. Ryan's world average of the lattice data for heavy 
light decay constants is in table |(| 

The lattice methods can be checked by computing the fo decay constant. 
The current experimental result for f Ds = 280 (40) (19) MeV |§. The CLEO-c 
experiment plans to reduce the experimental errors on fjj s to the few percent 



Focusing of a Parallel Beam 19 



Group 


method 


9* 


UKQCD [IK 1 
Abada et al. |lll] 


Lattice QCD 
Lattice QCD 


0.42(4)(8) 
0.69(18) 


CLEO 112| 


Experiment D * Dtt 


0.59 ±0.01 ±0.07 



Table 7. Summary of some results for g n 



level [ 1 07 1 to test lattice QCD. The actual comparison between theory and ex- 
periment will be ratios of matrix elements for leptonic and semi-leptonic decays, 
so that the test is independent of CKM matrix elements. 



In her review article Ryan |5j quoted the errors on £ from lattice QCD as 

-7 
-o 



£=1.15(6)±J (47) 

The first error in equation ^ is the statistical and systematic errors from 
quenched QCD. The asymmetric errors are from unquenching. It is instructive to 
compare the errors on £ with the experimental errors on Airid — 0.503 ±0.006ps. 
Although Am s has not yet been measured, it is expected to be measured to a 
few percent accuracy at the Tevatron. The final errors on | ^ | 2 will be limited 
by the theoretical errors on £. 

Unfortunately, during the last year the errors on £ have gone up again, based 



on some observations by the JLQCD collaboration ]108| |. Kronfeld and Ryan [109 
have suggested that a more realistic value of £ is 1.32 ± 0.10 rather than the 
estimate in equation |^]. 

The key problem is that the light quarks in the current unquenched lattice 
QCD calculations are not so light. Lattice QCD calculations are typically done 
at unphysically large mass parameters. Physical results are obtained by extrap- 
olating the results using effective field theories. 

The effective field theory for heavy-light systems contains the light particles: 
7T, K, and rj, and a pseudoscalar and vector heavy-light state The Lagrangian 
is written so that it is invariant under heavy quark symmetry and SU(3)l x 
SU (3)i? symmetry. This Lagrangian is for static quarks. The Lagrangian can be 
used to calculate masses of hadrons and decay constants in terms of the couplings 
in the Lagrangian. 

The most important coupling at tree level is the g w coupling that describes 
the D* — ► D + tt decay (suitably extrapolated to the heavy quark limit). Table [7] 
contains some estimates of from experiment and lattice QCD. 

The first loop correction to the decay constant has the form 

^TEf B = <2>(1 + Af q ) (48) 

where <P is the quantity with zero light quark mass and Af q represents the 
deviation from the chiral limit due to the finite size of the light quark mass. 

The problems with the chiral extrapolations of £ are due to the ratio of the 
decay constants, so consider: 

^-1=^-1 (49) 

JB 
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= - Sf d (50) 
= (m 2 K -ml)f 2 (n)+C (51) 

where f%{lj) is a low energy constant ol the effective field theory. The form of C 
is 

„ l + 3qi. 1 2 , 1 o, < m %s 3 2, 

° = l4# ( 2™^ ln( "# } + 4 TO ' ln( # } _ 4 m ' ln( 7# } (52) 

The equivalent expression for the bag parameters B has the coefficient 1 — 3g 2 

_ 2 

in front of the chiral logs, from the values of g v in table |7], the ln(^f-) term 
has a negligible effect. 

Until recently most lattice QCD calculations extrapolated £ with the C func- 
tion set to zero. For example, the MILC collaboration used linear and quadratic 
chiral extrapolations into their fits for their original results JIt). The JLQCD 
collaboration tried to fit equation [5l| to their unquenched data. 



Kronfeld and Ryan [ 109 1 noted that once the g„ in known, then the chiral 
log term in C is known. Hence, they used the lattice data that is essentially 
consistent with linear quark mass dependence and add the log term by hand. 
The problem with this type of approach is that it assumes that the current lattice 
data is in the regime where there are no higher order corrections to equation ^l|. 
The definitive answer for the value of £ will come from unquenched calculations 
with light quarks that explicitly see the chiral logs in fg. 

This "case study" demonstrates the importance of the parameters of the 
dynamical quarks to the computation of heavy-light matrix elements, particu- 
larly the masses of the sea quarks. This study also demonstrates that the use 
of quenched QCD to compute heavy-light matrix elements is coming to an end. 
The chiral structure of matrix elements in quenched QCD can be very different 
to that in unquenched QCD. 



5.1 Computation of form factors from lattice QCD 

One of the best ways to extract the | V c b | CKM matrix element from experi- 
ment is to use the B — > D*lvi semi-leptonic decays The differential decay 
, based on HQET ^ is 

^^^^^^^ (53) 

where K(w) is a known phase space factor and T(w) 2 a form factor. The value 
of w is the dot product of the velocities of the two heavy-light mesons. The 
expression in equation is based on heavy quark effective field theory 

As the masses of the b and c quarks go to infinity the normalization point of 
the Isgur-Wise function is obtained .F(l) = 1. The form factor at zero recoil is 
broken into the following 



■^(1) = VqedVa{1 + S 1/m 2 + ..) 



(54) 
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where t\qed is a perturbative QED factor and r\A is the perturbative matching 
factor between QCD and HQET. The &i/m 2 Q term represents the breaking of 
heavy quark effective field theory. A term of the form -=- is forbidden by Luke'e 



theorem [114]. 
The value of 



V c b I in the particle data tabic ]113| from B — ► D*lv decays is: 



cb \exclusive~ 



(42.1 ± l.l expt ± 1.9 t/leon/ )xl0 



-3 



(55) 



The theoretical error is dominated by the theoretical uncertainty in $\/ m 2 Q ■ I n 
the past <5i/ m ^ has been computed using sum rules and quark models. Without a 

systematic way of improving the results, 5% will be a lower limit on the accuracy 
of this CKM matrix element. 

In figure ^ I show a space-time diagram that is used to calculat e form factor s 



for semi-leptonic using the path integral. The initial lattice studies [115,116 117,118) 
mapped out the dependence of the B — > D form factors on w. The experi- 
mental dat a is taken at nonzero recoil w <> 1. However, the extrapolation 
to w = 1 [HE] is either done using a simple ansatz or using the results of a 



dispersion relation [ 12C , 121 



The Fermilab group [122] have concentrated on estimating ^i/ m ^ i n equa- 
tion ^4| They compute the matrix element in figure | with all the mesons at rest. 
By taking clever combinations of matrix elements they can get a precise estimate 
of the form factor at zero recoil. The mass dependence of the form factor can 
be mapped out by varying the masses of the heavy quarks. The matching to the 
continuum is quite involved. The final results [122] have errors with the same 
order of magnitude as other approaches. Their error includes a 10% estimate for 
unquenching. 

There are other ways that lattice QCD can contribute to the extraction 
of the | V c b | CKM matrix element from experimental data . There have been 
calculations of the semi-leptonic decays of the Ab baryon [123] fro m lattic e QCD. 
There have been two calculation of the mass of the B c meson [124 1231. Jones 



and Woloshyn [124] computed the decay constant for the leptonic decay of B c 
to be 420(13) MeV (the error is statistical only), using NRQCD at a lattice 
spacing of 0.163(3) fm. Whether these additional channels can help reduce the 
theoretical uncertainty of | V c b | to below the 5% level is not clear. 



6 Nonleptonic decays 

One of the main goals of the B physics experimental program is to check the 
CKM matrix formalism by measuring the CKM matrix elements many different 
[ 126] ,127|. For example, the experimental measurements for B + — > Tr + 7r° 



ways 



could be used to extract sin(2a) if the hadronic uncertainties could be controlled. 

The path integral in equation |l| is calculated in Euclidean space to regulate 
the oscillations in Minkowski space. This means that the amplitudes extracted 
from lattice QCD calculations are always real. Recently, there has been some 
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theoretical work (see [ 128 1 for a review) on the non-leptonic decays of the kaons, 
motivated by the attempts to compute the hadronic matrix elements for e'/e. 

There were some early attempts to study the decays D — » Kir on the lat- 
tice 129 13C , 131 1. These type of lattice calculations stopped when the theoretical 
problems with making contact with experiment became apparent. In this sec- 
tion I briefly describe some of the old work on the D — > Kir decays and provide 
pointers to the new theoretical developments. 

The correlator required is 

G(t) = (0 | (irK)(t)H eff (0)D(t K ) | 0) (56) 

where D(tx) is the interpolating field for the D meson at time tx and irK is 
the interpolating operator for the pion and kaon. The effective Hamiltonian is 

H eff = c + (/i)0; ont + c_{n)O cmt (57) 

where c±(/z) are perturbative coefficients. 

The diagrams for the Wick contraction of equation |5(] are in figure 0. Al- 
though the diagrams are more complicated to compute than those for leptonic 
or semileptonic decays they can be calculated on a supercomputer. The problems 
occur trying to extract the pertinent amplitudes from G(t) in equation |5^. 



As pointed out by Michael [132] and, Maiani and Testa |133| there is a 
complication with creating a pion and kaon state with definitive momentum. 
The operator 

7rK {t)=ir(p,t)K(-p,t) (59) 

has the same quantum numbers as irK with all possible momentum values. The 
ground state of the operator in equation ^9] will be the pion and kaon at rest. 
Hence, in the analogue of equation ^ for this graph, the required amplitude will 
not be the ground state. It is not easy to fit a multi-exponential model to data, 



although it is possible with a basis of interpolating operators [ 128 1 . 

Theoretical work, in the context of K — > irir decays has shown how to get 
matrix elements in infinite volume from matrix elements computed in finite 



volumes [134 ,135 1 . There have also been proposals 1 132 136] to introduce some 
model independence to extract the complex phases of the matrix elements. 

The methodology for non-leptonic decays will be further developed and tested 
on K — > irir decays, before any attempts are made at the decays D —* Kir, 



7 Conclusions 

The consumers of lattice QCD results need the error bars on current matrix 
elements to be reduced below 5%. The hardest error to reduce is from quenching. 
Improved staggered quarks look like they will be the first to explore unquenched 
QCD with light sea quarks. This should motivate the champions of other light 
quark formalisms to speed up their unquenched calculations. The techniques 
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that will reduce the error bars of heavy-light matrix elements will lie outside 
the domain of heavy quarks, in areas such as algorithms, improved computer 
hardware, and better grant writing. 

My own, admittedly biased view, is that quenched QCD calculations are 
now of limited use for lattice QCD calculations with heavy quarks. As every 
experimentalist I have ever met has held this view, I am sure it will prevail. 

The computation of matrix elements for two body hadronic decays still looks 
quite hard. Interesting things seem to be happening for kaon decays and in 
theory. It is not clear, whether these developments will be useful for non-leptonic 
decays of the B meson. It would obviously be a m ajor breakthrough if this 
problem could be solved, however Mark Wise's wise [ 127 1 words, on the career 
ending nature of working on nonlcptonic decays should be heeded. 

The computation of the QCD matrix elements for the heavy flavour program 
is a well defined task. If we can't compute them reliably, then we will have 
failed. We will have to admit that we can't compute anything from QCD outside 
perturbation theory from first principles. I hope this doesn't happen. 
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Fig. 2. Effective mass plot 
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Fig. 3. Spectrum of light hadrons from CP-PACS 
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Fig. 4. Spectrum of charmonium from CP-PACS 
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Fig. 7. Wick contractions for equation 



